Superfluid phases of quark matter, II: Phenomenology and sum rules 
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We derive sum rules for a uniform, isotropic superfluid quark-gluon plasma with massless quarks, 
first laying out the phenomenological equations obeyed by a color superconductor in terms of macro- 
scopic observables such as the superfluid mass and baryon densities, and the electric and magnetic 
gluon masses, and then expressing these quantities in terms of equilibrium correlation functions. 
From the transverse part of the long wavelength baryon current-momentum correlation function 
we derive an exact expression for the superfluid baryon density, and from the longitudinal part, 
an /-sum rule. From the transverse part of the long wavelength color current-current correlation 
function we derive the superfluid Meissner mass, and from the longitudinal part, the Debye mass. 
These masses constrain integrals of the transverse and longitudinal parts of the gluon propagator 
over frequencies, and provide self-consistent conditions for a solution to the gap equation beyond 
O i 1 weak coupling. 
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Dense degenerate quark matter is expected to exhibit color superconductivity in the color-antitriplet channel; the 
predictions are based on a weak coupling analysis of the gap equation (l]^ , as well as Ginzburg-Landau theory . 
£f) • In possible physical realizations in neutron stars and ultrarelativistic heavy- ion collisions, such a superconductor 
>- would be in a strongly coupled, color-singlet state. The equilibrium properties of such superconducting matter have 
ON i yet to be clarified in detail. In the strong coupling regime, Ginzburg-Landau theory delineates the possible phase 
diagrams near the critical temperature T c , but quantitative predictions for the fundamental parameters of the theory 
are lacking. 

In this paper we derive exact sum rules obeyed by the transverse and longitudinal momentum-momentum and baryon 
current-momentum correlation functions and the gluon propagator, D, in color superconductors. The sum rules, which 
are related to the linear response of the equilibrated many-body system to a current-inducing perturbation, connect 
the long wavelength behavior of the correlations with macroscopic observables 10 HJ|. Such observables include the 
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superfluid mass density, the superfluid baryon density, and the magnetic mass or inverse penetration depth, in addition 
to quantities such as the charge conductivity and the Debye screening length that play a role in the normal state. 
These sum rules act as self-consistency conditions that must be satisfied by approximate theories of thermodynamics 
and correlations in the superconducting state. To derive the sum rules we first set up the phenomenological equations 
obeyed in a relativistic superfluid plasma in terms of the macroscopic observables, and then turn to the expressions 
for these observables in terms of correlation functions. We consider a uniform, isotropic color superconductor of 
three-flavor massless quarks at finite temperature, T, and baryon chemical potential, fib, and use units h — c = 1. 



I. PHENOMENOLOGY OF RELATIVISTIC SUPERFLUIDS 

Let us first review the phenomenological equations obeyed by the baryon current and momentum density of a 
relativistic superfluid plasma in non-dissipative hydrodynamics, linearized about equilibrium with small velocities. 
Through these equations we derive the relativistic relation of the superfluid mass density, p s , and the superfluid 
baryon density, n s , quantities familiar in the nonrelativistic context fll2| . Consistency of the hydrodynamic equations 
dictates as well the form of the superfluid acceleration equation in a relativistic superfluid. The extension of the 
relativistic superfluid hydrodynamic equations, (^])-(||). (|l0|), and ( p^ ) below, to arbitrary velocities may be found in 

& 

The momentum density, g {gi = — Toj, the off-diagonal components of the stress tensor), is given in terms of the 
(small) velocities v s of the superfluid and v„ of the normal components by 

g = p s v s + p n v n , (1) 
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where p s is the superfluid mass density, p n is the normal mass density, and p s + p n = p + P, where p is the total 
mass density in the system at rest (the internal energy) and P is the pressure. That the superfluid velocity is an 
independent thermodynamic degree of freedom is a fundamental property of the paired state. In non-relativistic 
superfluids, p + P reduces to ran, where m is the rest mass of the carriers, of density n; relativistically one must 
retain the contribution of the pressure in the momentum density. Similarly the baryon current is given in terms of 
the normal and superfluid velocities by 

jb = n s v s + n„v„, (2) 

where n s is the superfluid baryon density, n n the normal baryon density, and n s + n n = rib, the total baryon density. 
The superfluid mass density and superfluid baryon density are closely related. As we derive below, 

p s =p b n s , (3) 

while the normal density obeys 

Pn = Mb«« + Ts, (4) 

where s is the entropy density. The latter follows from Eq. (||) together with the relation for the thermodynamic 
internal energy density, p = pb^b + Ts — P. Thus 

g = Pbjb + Tsv n . (5) 

The basic hydrodynamic equations for the superfluid are the equations of momentum and baryon conservation, and 
of entropy flow. In linearized hydrodynamics, g is driven by pressure gradients according to 

^+VP = 0. (6) 
at 



Baryon conservation reads as usual, 



£ + V-W (7, 



Since the entropy in a superfluid system is carried only by the normal fluid, in the absence of dissipation, the entropy 
density obeys 

^+V-(v„ s ) = 0. (8) 
Furthermore, to second order in the flow velocities, the conserved energy density, E — Too, is given by 

j-, Ps 2 i P n 2 m\ 

e = y v s + Y Vn + Pi ( 9 ) 

and the equation for conservation of energy is 

dE 

— +V-g = 0. (10) 

To derive Eqs. (^) and (Q), we explicitly calculate the time derivative of (^J), keeping only terms of second order, and 
use the above equations, together with the usual first variation, dp = pbdrib + Tds, and the Gibbs-Duhem relation, 
VP = n b Vp b + sVT, to find 

-T£- + V • (pbh +Tsv n ) = (j b - n b v n ) ■ Wp b + Ps(v s - v„) • (11) 

Identifying the energy current with the momentum density, we see that g is given by Eq. (|^), from which Eqs. (||) 
and (H) follow. 

In addition, the right side of Eq. ([ll]) must vanish identically. Eliminating jb there by means of Eq. (||) we find as 
a necessary condition for this term to vanish that v s obeys the superfluid acceleration equation, 

M6^+V^ 6 = 0. (12) 
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The superfluid acceleration equation follows directly from the fact that the baryon chemical potential and superfluid 
velocity are given in terms of the phase <f> of the order parameter by 



2 d( t> ™ 

r b = -m' (13) 



and 
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H b v s = V0. (14) 



The factor 2/3 is the baryon number per pair. We recall that the order parameter takes the form ^ a bfh{x) = 
(ip a f(x)ipbh( x )) = \^ abfh{x)\e l ^ x \ where ipaf is the spinor for quarks of color a and flavor /, and tpaf = ^S/ * s ^ ne 
charge-conjugate spinor in the Pauli-Dirac representation. In a non-relativistic system, the first /ib in Eq. ( |l2[ ) and in 
Eq. ( |l4| ) become simply the rest mass, to, of the carriers. 

An important consequence of Eq. ( |l4| ) is that the circulation is quantized according to 

2 

d£ ■ gA*b v s = 27w, (15) 
where the integral is around any closed path and v is an integer. 

II. MOMENTUM AND BARYON CURRENT CORRELATION FUNCTIONS 

We now derive the phenomenological superfluid densities p s and n s microscopically in terms of momentum and 
baryon current correlation functions, which characterize the linear response of the system to external disturbances, 
in particular here, a Galilean transformation. Consider the situation, following Ref. pl| , in which an infinitely long 
cylinder containing a color superconductor moves very slowly with uniform velocity v along its axis, which we take 
to be along z. We assume that the normal component is in equilibrium with the walls, so that v becomes the normal 
velocity, and that the superfluid component remains at rest. 

The response of the system to this motion of the walls can be described in terms of the transverse response to a 
static long wavelength perturbation, J c? 3 rg(r) • v, where the momentum density operator is 

9i = E ^l/HVAfc - f K b A m )^ b f + E( E ° x BQ )ii ( 16 ) 

af a 

the A" are the color gauge fields, with field tensors F£ v = d^A" — d v A" — gfa^A^AI as well as field strengths 
Ef = and Bf — —heijkFa > 9 ^ s the color coupling constant, and the A" b are the Gell-Mann matrices. 
The induced baryon current is given by 

(j b ) v = limy^ 1 (k ) 0)v, (17) 

k— >0 

where the baryon current operator is = -| Y] a f 'tpaf'ypafi and Xt^O*-'®) i s the transverse component of the baryon 
current-momentum density correlation function fll4| , 

xg'](M)= r^ flfoiMfrg) . as) 

■> J_ 00 2tt z-uj 

We write here, for general operators a(r, t) and 6(r, t), 

([o ) &])(k,w) = -i / d 3 (r-r') r d(t-t0e-*-( r - r ')e'^*-*')<[a(r,t) ) 6(r',t')]), (19) 



where (• • •) is the ensemble average at given T and /if,. The retarded commutator is given by taking the limit of z 
approaching the real axis from above in Eq. (|l^) . 

Comparing Eqs. ( p"7j ) and (|2|) we see then that the normal baryon density is given in terms of the transverse baryon 
current-momentum density correlation function by 
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n„= Hm X f 9l (k,0). (20) 

k— >0 

This equation is effectively a sum rule obeyed by Xt^ * Similarly, the induced momentum density is given in terms of 
the transverse momentum density-momentum density correlation function by 



\gg] 

where for complex frequency, z, 



lim X ^(k,0)v, (21) 

k— »0 



(22) 



Thus, Eq. (^l| ) with (Q) yields the transverse sum rule, 



limxJf 1 (k,0) = f M b n n + Ts. (23) 

k^O 



To derive the longitudinal versions of the sum rules (g0|) and ( p3| ) we suppose instead that the cylinder is still long 
but finite with closed ends. Then the superfluid component flows together with the normal component, leading to 
(jb)v = n b v instead of n n v. In this situation, the linear response analysis yields [jLll 



We thus obtain the /-sum rule 



<j h )v= limx^(k,0)v. (24) 

k^O 



X^ 9] (k,0) = n fc , (25) 



as |k| — > 0. Note that the /-sum rule ( p5|) can be directly derived, for general k, from the baryon conservation law 
(Q) and the equal time commutation relation, ([j b0 (r, t), g(r', t)]) = — mbV£(r — r'), where j b o — \ J2 a f Vaf^af' 
sum rule can be rewritten in terms of x^ 9 ' as 

lim x l l 9] (k, 0) = n b n b +Ts = p + P. (26) 

k^O 

To derive this result we calculate g from Eq. ([5J) to first order in v = v„ = v s , using (jb) v — n bV and the fact that 
the entropy term is explicitly first order in v. 

The four sum rules (po|), (|25|) , (E^), and (26) relate the total and superfluid baryon densities n b and n s , and the long 

wavelength behaviors of the correlation functions Xt 9 ^ Xl 9 ' j Xt^ > an d Xl 9 ' • At T = 0, as in ordinary superconductors, 
longitudinal first sound modes, the only low-lying excitations for colors and flavors involved in the pairing, contribute 

only to Xl 9 ' and Xl 9 ': an( i hence n s > jn|. 



III. COLOR PHENOMENOLOGY AND COLOR CURRENT CORRELATION FUNCTIONS 

Color superconductors have the property of screening out color magnetic fields, the color Meissner effect. Following 
the line of argument of Ref . JllJ] , we consider the linear response of the system to an applied static long wavelength 
transverse color magnetic field, Ag Xt (r) = A^ xt (k)e 4k ' r , where k-A^ xt (k) = 0. The external field Aj xt produces 
currents of various colors, /3, which in turn induce color fields A^ ld (r); the total color field is A* 3 = A^ xt + A^,. In 
the static long wavelength limit, the induced transverse color currents are given in terms of the total transverse color 
field by the London equation, 

<r T (r)) ALt = -(m 2 M UAP(r), (27) 

where Tom is the magnetic mass matrix, non-zero in the paired state. The inverses of its eigenvalues are the length 
scales on which color magnetic fields are screened in the superconductor. 

To linear order in A cxt , the long wavelength induced currents are given microscopically by 

Cf T (r))A ?xt = - lim x£ 7 (k,0)A^ xt (r), (28) 
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where \t is the transverse part of the color current-current correlation function, 



2tt 



and the color current operator for gluonic index a is 



abf 



The linearized field equation for A? nd , 



implies that the total color field is 



V x (V x Af ,) = (f T (r)), 



A^(k) = [(^)- 1 ]^A: xt (k), 
where eTg(k) is the static transverse color dielectric function, defined by 



[(s T )- 1 U(k) = 5 a p-x¥ (k,0)/|k| 



(29) 
(30) 

(31) 
(32) 

(33) 



One can write e^ ;3 (k) in terms of the screened correlation function, Xt (k, 0) = x^ 7 (k, O)e^g(k) (the irreducible 
bubble), as 

el /3 (k)=^ + x? /3 (k,0)/|k| 2 . 
The magnetic mass matrix is thus given in terms of \t by 



(34) 



(m 2 M ) af3 = lim £t ( k >°)- 



(35) 



The transverse screening lengths are closely related to the superfluid baryon density n s . We may see the explicit 
relation near T c , where we derived the equilibrium properties utilizing general Ginzburg-Landau theory |^,|l5|. Within 
color and flavor antisymmetric pairing channels having zero total angular momentum, even parity, and aligned chirality, 
which include the two-flavor and color-flavor locked condensates as optimal states, we obtain the London equation for 
the induced current densities (|2^) as (see Appendix and Jl5| ) 



<j QT (r)) A2xt = -K T (|) 2 Re {<Tr [((X a )*<f> + + <M Q ) ((^T<t>+ + F } A*(r) 



with coefficient 



K 7 



9n s 



(36) 



(37) 



Here (4>+)abfh is the pairing gap of the quark of color a and flavor / with that of color b and flavor h, and the subscript 
F denotes the gap calculated for the paired quarks lying on the Fermi surfaces; for a color neutral system, the Fermi 
energies reduce to a single value /i(,/3 as T approaches T c . The relation (37) is basically that obtained by Josephson 
p6[ , p s = A±(m/h) 2 \^\ 2 , for superfluid He II. Equations ( p7j ) and ( fj6| ) then yield the relation between the magnetic 
mass, the superfluid baryon density, and the order parameter near T c : 



l 2 \ 9 3 n s 



((A Q )*</>+ + </>+A Q ) ((A^)*0+ + 0+A' 3 ) t 



!6Aifc Tr((t)\(j) + ) F 
For the color-flavor locked phase, where (<ft+) a bfh = l^A^af^bh — &ah5bf), we obtain 

K T = 



3n s 



16/Xb|«A||' 



(38) 



(39) 
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and 



(m 2 M ) aa = (40) 



8fi b 

Similarly, for the two-flavor channel, where (<p+) a bfh = £fhs£abcd c (s, the strange flavor) 



16/X(,|d|j. 



and 



0, a = 1,2,3, 

{m 2 M ) aa = { 9g 2 n s /16fi b , a = 4, 5, 6, 7, (42) 
3g 2 n s /4^ b , a = 8. 

In contrast to the behavior of transverse color fields, color longitudinal fields are screened in both the normal and 
superconducting states. Longitudinal color charge correlations act to expel low frequency longitudinal color fields as 
in the response of nonrelativistic electron systems to an external longitudinal electromagnetic field JIo|]. We can see 
this behavior by simply replacing the transverse external field A^ xt above with a slowly varying longitudinal one, 
satisfying k x A^ xt (k) = 0, with time dependence e - l (^+ l v)t (here, uj ~ and r\ is a positive infinitesimal). Then, as 
in the derivation of Eq. (|32|), we obtain, for the total color longitudinal fields in a gauge where the scalar fields Aq 
vanish, 

A^k,u J ) = [(s L )-% 1 A2 xt (k,u J ), (43) 

where £^g(k, uj + irj) = S a p — Xjf '(k, uj + irj)/\k\ 2 , with the screened correlation function x^(k, uj + iff) = x^ 7 (k, u> + 
ii7)£^g(k, uj + irj), is the longitudinal color dielectric function. [Here it is more convenient to write x'rf = Xoo > m 
constrast to the definition of the longitudinal part in note Hj.] 

In the static long wavelength limit, the longitudinal correlation function reduces to minus the square of the electric 
mass tensor, 

limxf (k,0) = -(m|)^; (44) 

then A Q (k, 0) — > |k| 2 (m^ 2 ) Q ^Af xt (k, 0), so that the total longitudinal field is screened in the long wavelength limit. 
In the opposite limit of spatially uniformity, with slow variation in time, 

lim Ww/|k|) 2 xf (k,w + «j) - {UJ 2 P U, (45) 

lo^O k— »0 

where uj^ is the plasma frequency matrix; in this limit A a (0, uj) — > — uj 2 {uj^ 2 ) a pA^. Kt (0, uj). 

Near T c , the square of the plasma frequency is given by the sum of normal and superconducting contribu- 
tions, {ujp) a p — 3 a poj 2 + (m 2 M ) a p (see Appendix); in weak coupling, uj 2 — g 2 ^1 / '18ir 2 + g 2 T 2 /2 |l7) and Kt — 
7((3)nb/16n 2 T 2 iib where £(3) = 1.202... is the Riemann zeta function. The deviation of uj^ from 5 a puj p 

below T c comes from the fact, as clarified by Rischke et al. j?],[l9) for the two-flavor channel, that the condensate 
changes the color dielectric properties in such a way that the color superconductor is transparent to the low energy 
gluons having color charge associated with the two colors carried by a Cooper pair. In weak coupling, m^f reduces 
above T c to the usual Debye mass \/3uj p 5 a f3 [ p"7| . Below T c , however, m^f generally deviates from the normal Debye 
mass, due to the modification of the color dielectric properties by the condensate |@,^9|. Just below T c , {m E ) a fi 
behaves as 3w 2 (5 Q/ 3 — Z{m M ) a p (see Appendix). 

The relations of the Meissner and Debye screening masses, tum and the, to the color current-current correlation 
functions can be cast in terms of sum rules obeyed by the gluon propagator, 

Df v {vtyt') = -z<T[^(r,*)^(r',t')]>- (46) 

Since this propagator dominates the infrared structure of the gap equation ||, the sum rule constraints on the 
propagator are important to take into account in constructing a self-consistent solution to the gap equation through 
inclusion of polarization effects of the superconducting medium. 
The transverse propagator Dt, with spectral representation, 



G 



is related to the transverse part, x^(k, z), of the irreducible current-current correlation function by 

(D- 1 )^, z) = 6 a0 (z 2 - |k| 2 ) - xf(k, z). (48) 
Thus taking z — and the limit of small |k|, we derive the transverse sum rule, 

l m l / -7T- = {m M Up- (49) 

Similarly, the longitudinal propagator, Dl = -Doo i n the radiation gauge, has the spectral representation, 



Df(k,z) = ^S aP + I (50) 



kP .Loo 2tt 



1 f°° duBffaw 



z — u> 



and is given in terms of the longitudinal part of the irreducible current-current correlation function by 

(Dl 1 ) a0 (k,z)=6 aP \k\ 2 -xf(k,z). (51) 
Again in the static long wavelength limit, we find the longitudinal sum rule: 

,. / 1 . f°° dLjBf(k,U>)\ , _ 2A 

te [w Safi - L t,^^ ) = < m * >«*■ (52) 

In the normal state the sum rules (|4^) and (^) reproduce the results analyzed by Pisarski and Rischke Q in terms 
of the spectral representation up to one-loop order. The Meissner masses vanish in the transverse sector; Landau 
diamagnetism leads only to |k| 2 corrections, and not a nonzero screening mass in the limit k — + 0. The right side 
of Eq. ([l9|) is replaced by 6 a p lim^o |k|~ 2 to leading order in <?, as in Ref. ||. The longitudinal sector contains the 
usual Debye screening, which is characterized by {mr E ) a p = 2>u>p5 a p to leading order in g. 

Rischke explicitly calculated the zero temperature, low frequency, long wavelength limit of the color current- 
current correlation function for the two-flavor and color-flavor locked condensates to leading order in g. Repeating 
his calculations for three flavors, we find the screening masses (with no sum over a), 

21-81n2 2 



(m E ) aa = 3(m M ) QQ = — uj p , a = 1, . . . , 8 (53) 



for color-flavor locking, and 



(m|) QQ =w 2 , {m 2 M ) aa = 0, a = 1,2,3 
(ra|)s8 = 3w 2 , (m 2 M ) ss = 2w 2 /9 



(m%) aa = 2w|, {m\j ) aa = w 2 /3, a = 4, 5, 6, 7 (54) 



for the two-flavor channel in which only quarks of color R and G, which couple to gluons of a = 1,2,3, undergo 
BCS pairing. (The choice of the two colors involved in the pairing is arbitrary under the constraint of overall color 
neutrality.) The screening masses depend on color charge only in the two flavor condensate, because there the order 
parameter is anisotropic in color space in contrast to that in the color-flavor locked condensate. How static transverse 
screening influences the pairing gap compared with Landau damping of color magnetic gluons |pof depends sensitively 
on the gluon energies and momenta controlling the pairing gap |Q . 

Weak coupling calculations ignoring the effects of the superconducting medium yield the logarithm of the gap 



ln(A/// b ) = -37r 2 /\/2. g -51n 5 + .... (55) 

The weak coupling solution is not self-consistent in the sense that it satisfies the sum rules (^) and ( |52| ) with the 
masses m E = \fZu) p 5 a (3 and = 0, corresponding to the thermodynamics of the weakly interacting normal gas 
||. As Rischke || showed to one loop order, the superconducting medium significantly modifies the gluon self-energy 
from the normal medium value only in the energy range |o)| ^ A, which is not sufficient to change the logarithm of 
the gap from Eq. (]55|) up to subleading order in g. On the other hand, a self-consistent solution that satisfies the sum 



7 



rules ([l9]) and ( |52| ) with masses given by Eqs. ( [53"| ) and ( |5~i| ) would include contributions in all orders. The sum rules 
thus provide a check on approximate theories for the pairing gap beyond weak coupling. 

In summary, we have derived the transverse and longitudinal sum rules for a color superconductor that is uniform 
and isotropic in ordinary space. In doing so we have brought out the relation of the long wavelength behaviors of the 
momentum-momentum, baryon current-momentum, and color current correlation functions to macroscopic quantities 
such as the superfluid density and the Meissner and Debye masses, as well as the relevance of the sum rules for the 
gluon propagator to the self-consistent solution to the gap equation. The sum rules hold for any number of flavors 
and quark masses. It is straightforward to extend the present analysis for color current-current correlations to the 
case in which the electric currents coexist with the color currents; in this situation the electric currents modify the 
supercurrents through a mixing of color and electric charge (see, e.g., Refs. jl5| , ^l| ). 

We thank Dirk Rischke for helpful comments. This work was supported in part by a Grant-in-Aid for Scientific 
Research provided by the Ministry of Education, Science, and Culture of Japan through Grant No. 10-03687, and in 
part by National Science Foundation Grant No. PHY98-00978 and PHY00-98353. 



APPENDIX: THE GINZBURG-LANDAU REGION 



In this Appendix, we summarize the derivation of the London equation (|36|) near T c , for color and flavor antisym- 
metric pairing channels having zero total angular momentum, even parity, and aligned chirality; full details will be 



given in Rcf. 



In the presence of applied weak color fields, A^(x), of very long wavelength and low frequency, the 



resultant gradient of the order parameter adds a small correction to the homogeneous part of the Ginzburg-Landau 
free energy derived in Ref. ||. Up to second order in the gap, this energy correction is 



Q g = i^ T Tr[(A^+) t D^+] F + ^K L Tr[(D <t>+yD (f>+] F) 

^0+ = - \ig[{\ a )* 4>+ + <j>+X a ]A<^ with the total color fields A°. 



(56) 



where the covariant derivative is D^(j) + = — ^ig[(X a )* (f)+ + 0+A Q ]^ with the total color fields A^. The time 

dependence of the gap is, by definition, measured with respect to that in the equilibrium phase, i.e., with respect to 
the phase factor e - 2l Mf>*/3 arising in the order parameter from Eq. ( |l3| ) (see Ref. ||). As we show here, the coefficient 
Kt is given by Eq. ((3?]). The coefficient Kt is not necessarily equal to Kt since Lorentz invariance is broken in a 
many particle system. In weak coupling, Kt reduces to 3Kt 0- As required, £l g is invariant under global U(l) gauge 
transformations and flavor rotations, as well as under local color SU(3) gauge transformations. 

To derive the dependence of the coefficient Kt on the superfluid baryon density n S: Eq. (|37]), we consider the 
situation in the absence of color fields, in which the pairs move uniformly with small constant velocity v s and the 
normal fluid remains at rest. The phase factor of the gap in the fixed frame transforms by (j>+ —> e lP ' r <f)+, where P is 
the total pair momentum. The total momentum of the superfluid is then 



3n s P/2 = p s v s 



(57) 



In this situation, 



n g = -^ T P 2 Tr(0+^) F . 



(58) 



Following de Gennes [£2|] and Wolfle |2^] , we then obtain the baryon current density from the usual canonical equation 
for the gradient energy density ( 58 ) as 



2 JSk = l K tTt{ ^ 

3 5P 3 w+s 



-)fP- 



(59) 



Using Eq. ( |57| ) to eliminate P, and Eqs. (^) and (||), we derive Eq. (37). Note that the extra time variation of the 
gap leads to terms in fl g of order P 4 , which play no role here. 

The color current densities induced by the applied weak color fields near T c are composed of the superfluid and 
normal contributions, and The superfluid color currents can be calculated as 



ls SA?. 



-[Kt + S^Kl - #r)]f Im {Tr [((A Q )>+ + 0+A Q ) t d^+\ J 

+ [K T + S^ (K L - K T )} (|) 2 ^Re{Tr [((A")>+ + <t>+\ a ) ((A^)*0+ + <M") + ] f } 



(60) 
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Equation (60) reduces to the London equation ( p6| ) for applied static transverse color fields when the spatial variation 
is of sufficiently long wavelength that we can ignore the term containing c^0+. 

We conclude this Appendix by considering the linear response to applied longitudinal color fields near T c . For 
static color fields of very long wavelength the induced color density jf° can be obtained from Eq. ( |60| ) as jf° — 
(Kl/ KT){m 2 M ) a i3A^ where m 2 M is given by Eq. (B8|). In weak coupling, this color density, together with the normal 
contribution — 3uj 2 A a0 , leads to the square of the Debye screening mass matrix, {m 2 E ) a p = 3aj 2 5 a/ 3 — 3(m^-) Q| g. 
For uniform color fields varying very slowly in time, we calculate the induced supercurrent from Eq. ( |60| ) as j™ — 
— {rn 2 M ) a i3 A' 3 , where m 2 M is again given by Eq. (|38|). Combining this supercurrent with the induced normal current 
—uj 2 A. a , we obtain the square of the plasma frequency matrix, (uj 2 ) a p = 8 a puj 2 + {rn 2 M ) a p. 



B.C. Barrois, Nucl. Phys. B129, 390 (1977). 

D. Bailin and A. Love, Phys. Rep. 107, 325 (1984). 

D.T. Son, Phys. Rev. D 59, 094019 (1999). 

D.K. Hong, Phys. Lett. B 473, 118 (2000); Nucl. Phys. B582, 451 (2000). 
T. Schafer and F. Wilczek, Phys. Rev. D 60, 114033 (1999). 
R.D. Pisarski and D.H. Rischke, Phys. Rev. D 61, 074017 (2000). 
D.H. Rischke, Phys. Rev. D 62 , 034007 (2000); 62, 054017 (2000). 



D.H. Rischke, tiucl-th/0103050 



K. Iida and G. Baym, Phys. Rev. D 63, 074018 (2001). 

D. Pines and P. Nozieres, The Theory of Quantum Liquids, V. I (Addison- Wesley, Reading, 1989). 

G. Baym, in Mathematical Methods in Solid State and Super-fluid Theory, R.C. Clark and G.H. Derrick, eds. (Oliver and 
Boyd, Edinburgh, 1969), p. 121. 

L.D. Landau and E.M. Lifshitz, Fluid Mechanics (Pergamon, Oxford, 1987). 

I.M. Khalatnikov and V.V. Lebedev, Phys. Lett. 91A, 70 (1982); V.V. Lebedev and I.M. Khalatnikov, Zh. Eksp. Teor. 
Fiz. 83, 1601 (1982) [Sov. Phys. JETP 56, 923 (1982)]. 

Fhe decomposition of a tensor Xii (k, z) into its longitudinal and transverse parts is defined as usual by Xij = (&j — 
kikj /|k| 2 )xT + (kikj /\]i\ 2 )xL- The transverse response is realized here by taking k orthogonal to z, e.g., along x. 
K. Iida and G. Baym (to be published). 
B.D. Josephson, Phys. Lett. 21, 608 (1966). 

G. Baym and SA. Chin, Nucl. Phys. A262, 527 (1976); M.B. Kislinger and PD. Morley, Phys. Rev. D 13, 2765 (1976); 
D.J. Gross, R.D. Pisarski, a nd L.G. Yaffe, Re v. Mod. Phys. 53, 43 (1981). 



I. Giannakis and H.-C. Ren, tiep-ph/0108256 



D. H. Rischke, D.T. Son, and MA. Stephanov, Phys. Rev. Lett. 87, 062001 (2001). 

G. Baym, H. Monien, C.J. Pethick, and D.G. Ravenhall, Phys. Rev. Lett. 64, 1867 (1990). 

E. V. Gorbar, Phys. Rev. D 62, 014007 (2000). 
P.G. de Gennes, Phys. Lett. 44A, 271 (1973). 
P. Woffle, Phys. Lett. 47A, 224 (1974). 



9 



